ON THE BOLTZMANN DIFFUSION OF TWO GASES 



ALEXANDER SOTIROV AND SHIH-HSIEN YU 

Dedicated to Professor Tai-Ping Liu on the occasion of his 60th birthday 



Abstract. We study the Boltzmann equation for a mixture of two gases in one space dimension with 
initial condition of one gas near vacuum and the other near a Maxwellian equilibrium state. A qualitative- 
quantitative mathematical analysis is developed to study this mass diffusion problem based on the Green's 
function of the Boltzmann equation for the single species hard sphere collision model in 1131 . The cross 
species resonance of the mass diffusion and the diffusion-sound wave is investigated. An exponentially 
sharp global solution is obtained. 

1. Introduction 

The Boltzmann equation for a mixture of two gases in one space dimension takes the form of the following 
system: 

'd t F A + ed x F A = Q AA (F A ,F A ) + Q AB (F A , F B ), 

< 

d t F B +ed x F B = Q BB {F B ,F B ) + Q BA {F Bl F A ). 
The right hand side consists of the usual collision terms which for X,Y £ {A, B} are given by, [5]: 



Q XY (F x ,F Y ) = a XY [ [ (f* x F Y ' m -F x F Y *M-£< 

JS+ JR 3 



) • n\d^dn 



where F Y ^ = F Y (x,t, F Y * = F Y (x,t,^), and F' x — F x (x,t,£'), with S the unit sphere and S + are 
n such that (£ — £*) ■ n > 0. The constant factor a XY denotes the sum of the radii of molecule X and 
molecule Y. The post-collision velocities £' and are given by the formulas: 

r = €+ ^ (fo-fl-n)n, 
m x + m Y 

m x + m Y 

where m A and m B are the molecular masses of the two gases. These formulas are obtained from the 
conservation of momentum and energy at collision and the assumption that forces only act normally to the 
surface. 

The studies on gas mixture in terms of Boltzmann equations were initiated in |17j to study the condensation- 
vaporization problem for mixture of vapors of different species. Many interesting physical problems such 
as "ghost effect" and "Knudsen layer" for gas mixture have been investigated, [Tlll6|. 

This research is also interested in phenomena related to vapor-vapor mixtures. Its focus is on the basic 
phenomenon of mass diffusion of a finite amount of gas A dissolved into the surrounding gas B. From this 
consideration, we write 

{F A = fyiMj| for gas with finite total mass amount, 
F B = ffiM^ + M B for gas in the background. 

Here, and M B are the Maxwellian states 
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2RT J 
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m A 



/m B J y/{2TTRTf ~" r \ m B 2RT ) 
Here T is the temperature of the mixture, p B the density of gas B, R the gas constant, and we have 
assumed that the average velocity is zero. We also assume that 

RT = 1. 



The two Maxwellian states satisfy 

fQ^(IVU,IVU) = 0, Q bb (M b ,M b ) = 0, 
1 Q AB (M A , Mb) — 0, Q ba (M b ,Ma)=0. 



(1.1) 



A mixture of rarefied gases is important for many high technology industrial applications. A classical 
approach to the study of gas mixtures is through a Navier-Stokes model. However, to use the Navier-Stokes 
model one needs to obtain the diffusion coefficient through formal calculations of a transport equation 
which are carried out easily only for either small or large ratio of the molecular masses |11). This makes 
the foundation of the Navier-Stokes approach uncertain, especially since the Navier-Stokes equation is 
valid for fluid near thermo-equilibrium state. Often in applications the state of the gases is not necessarily 
close to such a thermo-equilibrium. The validity of the Navier-Stokes approach is also questionable if 
there is some physical boundary encountered. Such issues on validity for Navier-Stokes in the application 
to rarefied gases had been extensively studied by an engineering school in Kyoto University initiated by 
Professor Sone, |T2]- Thus, it is desirable to study the problem entirely through a fundamental physical 
model based on a kinetic equation. 

The system for (f A ,f B ) is 



where 



(1.2) 



Lab^a= 



L AB f A + r AA (f A ,f A ) + r AB (f A ,f B ), 

Lf B + Y BB (fg , fg) + L BA f A + T BA (f B ,f A ), 



J B A 



Mb = 



fA = 




) + q bb (m b j b Vm b ) / 

xJyVm^) for X,Y e {A,B}. 



Here, Y XY are all quadratic nonlinear terms. 

When both the total mass of gas A and the perturbations in gas B are sufficiently small, in general the 
basic time asymptotic behavior of the nonlinear solution is governed by a linear equation (up to a large 
time scale). It is reasonable that scent diffusion in a gas is a fast time scale phenomenon independently of 
whether the density of scent is strong or weak. Therefore in great generality the following decoupled linear 
Boltzmann system 

(1.3) d t g + ed x g = L AB g, 

(1.4) d t h+ed x h = Lh + L BA g, 

will give a refined quantitative description of mass diffusion. It can be written in a matrix form as 

'eo x - l ab o 

-L BA ^d x - L 

We will call 1|1.3[) the linear Boltzmann diffusion equation, and we will call the system (|1.5f) the linear 
Boltzmann system. There are several features which provide a qualitative difference between the problem 
at hand and the single gas flow or flow in a mixture where both gases are near a nonzero equilibrium state: 



(1.5) 



0. 
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First, a purely diffusive behavior is dominant in the gas near vacuum and so we have qualitatively 
recovered the behavior already expected from a classical fluid treatment. Acoustic waves are not present 
in the first gas, such waves however appear in the second gas. 

Second, as discussed in the paper [I] , the technique that would be natural in dealing with the problem 
of two gases, both near nonzero background states, would require to work in a functional space of pairs of 
functions (f^fs) with an inner product that depends on the densities pa, Pb- 

((fA,fs), (Ha, hs)> = Pa{U, M + Pb^b, hs>. 

This inner product will yield a self-adjoint operator on the space of pairs of functions on which a spectral 
analysis can be performed and a theory similar to that developed in ^3] can be constructed. Notice however 
that if one of the gases is near vacuum (so pa = 0) the inner product becomes degenerate and it becomes 
impossible to work with the operator on pairs of functions. This is the problem that is addressed by this 
work. 

Third, the characteristics for the macroscopic equations in the system <|1.3ll . <|1.4fl coincide along the 
direction of mass diffusion. A further physical property of the collision operator, namely a microscopic 
cancellation representing the conservation of the total mass is employed to resolve this resonance. 

To see the generic feature of mass diffusion the following form of the initial data (g, h)| t= o = (gj n , hj n ) 
will be sufficient: 



(1.6) 



gin(x,Q = h;„(a:,0 = for \x\ > 1, 
sup(l + |£|) max(|h in (a:, |h 

in («£ 5 

01) <i. 



\x\<l 

. £eR 3 



Theorem 1.1 (Main Theorem). Suppose the initial data are as in H1.6(l . Then the solution to (|1.3(l , ()1 .4|) 

satisfies the estimates 



\\h(x,t)\\ 



where C > 0, 



0(1) 
0(1) 



(x- \ 2 t) 
e C(l + t) 



-{t+\x\)/C 



e c<i+t) 



1 v^r+*j 



+ e 



-(t+\x\)/C 



10 for \x\>2\Xi\t, 
'\e~ t2/ /or Izj < 2|Ait| 



Ai — — J |, A 2 — 0, A 3 — ^ ., 



l|f(M)IU| = 



f(x,i,C) 2 ^ 



Within a finite Mach number region |x| < 2|Ai|(l + t), there is a tail decaying exponentially in time: 

0(l)e~ t5 / c . This is a purely kinetic phenomenon (such a tail should not appear in a linearized equation 
modeled properly after the Navier-Stokes equation). However, since this tail decays exponentially in time, 
it is difficult to recognize the difference in the solutions of kinetic equations and classic continuum fluid 
mechanics when the time scale is large. 

The basic principle used to develop the theory is the "separation of scales", which is a concept 
created in ^5]. I n a Long Wave-Short Wave decomposition and a Particle- Wave decomposition were 
introduced. These two decompositions are designed for solutions with different physical characteristics. For 
example, when the solution is in a region corresponding to the "small time scale" or the "high Mach number 
region" (large space scale), the behavior is particle-like. In this case the Particle- Wave decomposition is 
used to determine the dominant behavior. When the solution is in a region corresponding to large time 
scale or in the finite Mach number region, the behavior of the solution is rather close to continuum flow. In 
this case Long Wave-Short Wave decomposition is used to analyze the solution. These two decompositions 
can be matched to each other with the difference decaying exponentially fast in time in the norm |j • ||x, oo (z, 2 ) 
Note there is no exponentially decaying structure in the x variable. This is an instance of the uncertainty 
principle from harmonic analysis. The Long Wave-Short Wave decomposition allows one to work in Fourier 
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space and there use the microscopic cancellation to resolve the resonance problem due to the waves along 
the characteristic x = 0. 

In Section II we will describe the physical problem clearly and review some of the basic theory of the 
Boltzmann equation for a single gas relying on E| • in Sections III and IV we establish the spectral 
analysis for the Boltzmann Diffusion equation and continue to establish its Green's function. 

In Section V, a microscopic cancellation and a "scale separation" are discovered for the gas mixture 
problem. The main theorem is proved. 

We would like to note that the above result allows an extension to a solution of the fully nonlinear equa- 
tion provided the initial data is sufficiently close to equilibrium. The first author has carried out such an 
extension as a part of his dissertation in the case when the molecular mass satisfy the relation > ■ 
While the proof involves lengthy computations and a number of genuinely nontrivial difficulties had to be 
overcome, the main techniques are similar to those used for the nonlinear problem in |13j . 

Acknowledgment. The authors would like to thank Professor Tai-Ping Liu. His penetrating vision 
on both scientific and mathematical problems has been influencing the second author from the very begin- 
ning of his scientific career. Both authors have greatly benefited from his teaching, his sharing of scientific 
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sixtieth birthday, the first author being a current graduate student of his, and the second one of his former 
graduated students, a friend, and a long time collaborator. 



2. Preliminaries 
2.1. A Physics Problem from Rarefied Gases. 

The primary content of this problem is to study the diffusion of Gas A with finite amount of mass in Gas 
B which occupies the entire space of a one-dimensional pipe in terms of a system of Boltzmann equations. 
Since this is a linear coupled system, the linear mass diffusion of Gas A will have some resonance from 
the background particle diffusion of Gas B. The linear resonance may also have a significant effect to the 
acoustic waves propagating in Gas B. It would be interesting to determine which one of Figure 1 and Figure 
2 below gives a correct picture of the mass diffusion for rarefied gases (notice the difference in spatial decay 
in the two figures). 



Diffuse of Gas A in a Pipe full of Cas B j-o 



Diffuse of Gas A in a Pipe full of Gas E - 1 -° 



Sound Wave in Cas E 

Figure 1. 



Region confining 
Concentration of Gas A 



Sound Wave in Gas I 



Sound Wave in Cas B 

Figure 2. 



Region confining 
Concentration of Gas A 



Sound Wave in Gas B 



Here, Figure 1 is a picture based on common experience about gases close to an equilibrium state. 
However, our knowledge for a gas near thermo-equilibrium is not necessary valid for a rarefied gas. Figure 
2 is a picture adopted from the consideration of resonance and coupling of two different gases. One might 
conjecture that there is a tail algebraic in space after the sound wave. 

In general in physics, one may need to set up a physical experiment to judge which figure is close to 
the reality of a rarefied gas. Unfortunately, to setup an experiment for this problem could not be easy in 
practice. Since the question is about the coupling of mass diffusion to sound wave propagation in a rarefied 
gas, it would require a rarefied gas environment with a large length scale. 
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Another alternative in physics available nowadays is a numerical simulation from first principles which 
may give some information about the truth. However, the requirement for a larger space size and the 
number of particles 1 makes a reliable numerical simulation very difficult. 

Thus, a quantitative and qualitative mathematical analysis would be a good alternative for this problem. 
Such a mathematical analysis would benefit other physics problems in the field of rarefied gases. 

2.2. Mathematical Theories on Boltzmann equation. 

For the convenience of the reader, in this section we include the necessary theorems established for a 
linear Boltzmann equation: 

(2.1) d t b + ed x b = Lb, 

where L is the linear collision operator given in (|1.2|l . The linearized collision operator L is symmetric on 
the student Hilbert space L 2 given by 

f(h,g)= f h(0g(£Rforh,ge4 

[\M\q = VM- 
For a given function g(x, t, £), we have the following notation: 



Expression 


Definition 


Comment 


llglUi 




non-negative real-valued function in (x, t) 


l|g|Uj(i|) 


(ff g(x,t,0 2 didxV 

\JR JR 3 ) 


non-negative real- valued function in t 


l|g||iT*(i|) 


(f f ]Tfe(z,i,oi 2 ^) 

\Jr Jr 3 l=0 ) 


non-negative real- valued function in t 


llelk«(i|) 


SUp||g|| L 2 


non-negative real- valued function in t 


lllslll 


sup |g(a;,i,e)|(l + |el) J 

££R 3 


non-negative real- valued function in t 



For the linearized collision operator L linearized around M = Mb , the null space of L is a five-dimensional 
vector space with orthogonal basis Xi, i = 0,l,2,3, 4, 0J: 

ker(L) = span{xo,Xi,X2,X3,X4}, 
Xo = MVa 

^ = eM 1/2 for i= 1,2,3, 
X4 ^-l f (|e| 2 -3)M 1 /2. 

Since we are interested in planar wave propagations, we restrict the function space to: 

(2.2) i|^{gei oo (IR 3 )|(e 2 M 1 / 2 ,g)=0, (? 3 M 1 / 2 ,g) = 0, ||g|U|<oo}, 

Macro-Micro Decomposition 

This space is invariant under the full Boltzmann operator Q BB and the linearized operator L. We decom- 
pose the Hilbert space L| = ker(L) © ker(L)- L : for any g G L|, 

!g = P g + Pig(= go +gi), 
Pog = (xo,g)xo + (xi,g)Xi + (X4,g)X4, 
P lg = g-P g. 



The particle density of a rarefied gas at 10 6 atm, 273° K is in the order I0 16 /£. 
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The characteristic information of the Euler equations is connected to the operator Po^ 1 on PqL^, \V2\ : 

dim(PoLt) = 3, 

Po^Ei = XiEi for i = 1,2,3, 

{Ai = - A 2 = 0, A 3 - V5?3} , 



(2.3) 



Ei 
E 2 
E 3 



-,Xi + X4 , 



iXi + Xij , 
,(Ei, Ej) — Sj (Kronecker's delta function). 

c = I Ai| : Speed of sound at rest. 
For the hard sphere model we consider, the linearized collision operator L is of the following form, 9 : 



(2.4) 



Lg(0 = -HZMO + Kg(0, 

Kg(0 =L,tf(£,£.)g(£.)rI£. 



■ exp 



■ exp 



2e"4-+2 + 



(l£| 2 + l£*l 



> ^(1 + 



Lemma 2.1. There exists positive constant i/ such that, for any h G L|, 



(P h,P h) = ^(E J ,h) 2 , 

3 

(P h,e 1 Poh)-^A,(E J ,h) 2 , 
[(Pih,LPih)<-^(Pih,(l + |e|)Pih). 

This lemma follows from direct computations and Carleman's theory on the negative definiteness of the 
operator L on Range(Pi), 0. 



Estimates on Green's function. ^31 

The Green's function <E(ir, t, £, £*) for the initial value problem l|2.1(l constructed in ^3] plays a fundamental 
role for studying the Boltzmann diffusion equation. 

It has also been shown in ^3] that G(a:, also satisfies the backwards equation: 

(-0 T - £*d y - L)G(x -y,t- r,£,£*) = 0, 
G(x-y,0,Z,&=6 1 (y-x)6 3 

and the function 

b(£, £,£)=/ / G(x-y,t,£,£*)b(y,0,£*)d£*d?/ 
Jr Jr3 

solves the initial value problem Ij2.1|l . 

Without spelling out the parameters £ and in the Green's function G(x, t, £,£*), one can treat G(x, t) 
as an L?-operator- valued function in (x, t): 
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For given (x, t) e 



G{x, t) ; h e L\ i — ► G(x, f)h £ L|, 



G(M)KO 



Analytic Spectral Decomposition 

Take Fourier transformation of 12.111 in x variable, then the solution b(x,t) can be expressed 

b(fc,i) = e ( - 4fc « 1+i) *b(fc,0) for fcei 

Lemma 2.2 (Ellis-Pinsky). 

Assume satisfies the restriction l|2.2(l . 

(%) For any ko > if \k\ > kq > 0, k 6 M i/iera i/iere exists e(fco) > suc/i i/iai z/cr zs an eigenvalue of 
—ikt; 1 + L then Re{a) < — e. 

(ii) There exist kq and S such that for \k\ < kq the spectrum with \cr(k)\ < 5 consists of three points 
o~i(k),o~2(k),o~3(k) which converge to the origin as k tends zero and they , together with the corresponding 
eigenvectors are analytic functions of k for \k\ < kq. 

(Hi) More precisely, the expansions of the eigenvectors £j(k) and eigenvalues <Jj{k) are as follows: 



(2.5) 



(2.6) 

and, 



o-i (*) =-i\ 1 k + A 1 \k\ 2 + 0{l)\k\ 3 , 

<r 2 (k) = -i\ 2 k + A 2 \k\ 2 + 0{l)\k\ 3 , 

a 3 (k) = -iX 3 k + A 3 \k\ 2 + 0(l)\k\ 3 , 

A, =(P 1 f 1 E J ,L- l P 1 e 1 E i )<0 ) 

A 3 = Ai, 



{-i£}k + L)&j(k) = crj(k)ej(k), 
(ej(k),e k (k)) = 6 jk , 
^(k) = E l + ke' l (0) + O(l)\k\ 2 , 

3 



3 = 1 



(Xj - X k ) 



for j ^ k, 



4 = o, 

e£ =ii- 1 PiC 1 E fe . 



Note that the analyticity property of o~j(k) for k around follows from theorems in 

Lemma 2.3 (Liu-Yu). There exists k,q > such that for all \k\ < kq the operator e^ jfc ^ on L 2 - can 

be decomposed as 

(2.7) 



{-ikt}+L)t 



^V^e^fc) <g> (e,(fc)| +e(- lfe « 1+L ) t nr: 



so </ia< 
(2.8) 



e ( -*= 1+i) *n^ < Ce"* /C /or some C > 0. 



Here, the operator e.,(fc) (g> (e.,(fc)| is given as follows. 

For any j and k in L 2 , the operator j <g> (k| on L 2 is defined as follows: 

(k|g= (k,g)j. 

Now taking inverse Fourier transform in the generalized sense and writing as a convolution we get for 
the solution: 

g(x,i)=G*g m = -L / e lkx+ ^ lk+L ^g m (k)dk= f G{x - y> t)g in (y)d y> 



8 



ALEXANDER SOTIROV AND SHIH-HSIEN YU 



where 



G(x,t) 



s ifcx+(-ie 1 fe+L)t rffc 



and we accept the last identity as the definition of Green's function G(x, t). G(x, t) is a generalized operator 
valued function of x with values operators on L| = L 2 (R 3 (£)) since we obtained it by taking inverse Fourier 

transform of the regular operator valued function e' - ^ fe+i ) t which as a function of k does not necessarily 
have a convergent Fourier integral. By pointwisc estimates we mean estimates of its operator norm as a 
function of the points x and t: but since the function is generalized of course one needs to identify and 
subtract its singularity for such estimates to make sense. However at this stage we will only look for 
estimates once Green's function is convoluted with a bounded function with compact support in the x 
variable, so pointwise estimates make sense; i.e. we will look to estimate ||G<7j n || where gi n is as in l]1.6[l . 
Long Wave-Short Wave Decomposition. 



One can decompose the semi-group G(x,i) 



'2tt 



Jxk+i-ik^ 1 +L)t 



dk into 



(2.9) 



G L (x,t) 
G s {x,t) 



/2tt 
1 



J.xk+(-ik^+L)t^ 



|fc|<K 



|fe|>KQ 



ixk+{—ikt}-\-L)t 



dk. 



The Green's function is obtained in terms of the Fourier variable k without detailed spectral properties 
of the operator (— ik^ 1 +L) described when \k\ > Kq. Thus, it not possible to obtain ||G,g(a:, t)\\ L 2 pointwise 
through Fourier analysis alone. This is a matter of uncertainty principle. 

Particle- Wave Decomposition. [T5| 
Consider an initial value problem 

(d t b + ed x b = Lb, 

\ b(x,0) = g in , 



(2.10) 
where gi n satisfies 
(2.11) 




for \x\ > 1, 

gin||L~(L|° 3 ) = 



1 < oo. 



Lemma 2.4 (Liu-Yu). For each given j £ N, there exist functions P j (x , t) andWj(x,t) so that the solution 
b(x, t) of i|2.1U[) with the given initial data in (|2.11l) satisfies 

'b(x,t) =F j (x,t)+W j (x,t) for all x e R,t > 0, 
F j (x,0)=g in (x), 

\\Fj(x,t)\\ L i < Cje-^+V^J for allxeR,t> 0, 

ll^illff j (L 2 ) — for all t > (uniformly bounded in time), 
where Cj > is an universal constant. 

Remark 2.5. The existence of such a decomposition is due to both the introduction of an essential kinetic 
equation (4.2) in |13| to resolve the singularity in the initial data into source terms which are smooth in £ 
and to the discovery of Mixture Lemma in |13|. ■ 

Theorem 2.6 (Liu-Yu). For the L^- operator-valued function G(x,t), Lemma \2,.J\ is still valid, i.e. 
exist Fj(x, t) and Wj(x, t) such that 



there 



(2.12) 



G{x,t) = F j (x,t)+W j (x,t) for all (x,t) e R x R+, 

F,(x, t) = e-^5(x - et)5 3 (£ - e.) +h(x, t) + j 2 (z, t) + • 

||ji(ss,t)|| z a < Cje-lW+W ' for allxeR, « = 2, 3, • - • ,j,€ 



+ ij(x,i), 
,i>0, 



\(d t +ed x -L)Fj\\ Hm) 

I^jWhUl 2 ) — Q f or a ^ * — Q( un tf orm ly bounded in time), 
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We recall the following theorems in on the Green's function G for the initial value problem. 

Theorem 2.7 (Liu-Yu). For the bounded initial data gi n with compact support in x, H2.11JI . there exist 
positive constants Cq and C\ independent of gi n such that 

for\x\ < let, 



HGW^IU? = 0(1) 



E 



e C (t+1) 



1 y/fjk + Tj 



J*Pig in (x)|| L 2, ||PiG t g m (a:)|| i 2 



0(1) 



e c <t+i) 



-[ + 



E 



-(|x|+t)/C! 



-(|x|+t)/d 



(2.13) 



3 



\Px&P lgin {x)\\ Ll = 0(X) \\\g m \\\ 



~[ (t+1) 3 / 2 



-k | &in 



Mil 



O(l) 



e c (t+i) 



-(|x|+t)/d 



V¥TT) ^ 3 (t + i) 



^ | e -(l- ' ' ■ 



/ 



/or |x| > 2ct 



|G*gin(a;)||i= < O |||gi 



-(M+t)/d 



where 



(x) = / G(x - y,t)g in (y)dy. 



Remark. This theorem is stated as Theorem 5.6 in [2]. The first work to employ complex analysis to 
obtain such exponentially sharp estimates is in f° r the Green's function of compressible Navier-Stokes 
equation. [x] 

Theorem 2.8 (Liu-Yu). For the bounded, compact-supported initial data gi n , (|2.11|l . there exist positive 

constants C\ and Co such that for all x £ R the main fluid part satisfies 

(2.14) 



[Ej ® (Ej| G* E fe 



=fc| fen](a:)IU t -, = 0(l)|||s*. 



5: 



__(x-A i £)- 
g C (t + 1) 



3 



, v — 



e c <t+i) 



-(|x|+t)/d 



onrf t/ie non-fluid parts have higher rate of decay in time: 



|[G t Pig i „(a;)|| i| . 8 =0(l)|||g in ||| E-(I 



e c <t+i) 



|PiG*g in (x)|U ? =0(l)||| gi . 



|PiG t P lgin (a;)|| L | ?3 = 0(l)|||g i 



~t (! + *) 

Q (^-A^t) 2 

* e - C (t + 1) 

(x-Ajt) 2 

e o (t+i) 



(l + <) 3 / 2 



-(|x|+t)/Ci 



-(M+t)/Ci 



-(|x|+t)/d 



3. Spectral Properties of the Linear Boltzmann Diffusive Equation 
3.1. On the Cross Species Linear Collision Operator. 

Similarly to the linear single species collision operator, the operator L AB is symmetric in L|, however now 
dim(ker(L AB j) = 1. One has the following lemma: 
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Lemma 3.1. 

i) The linear cross species collision operator L AB is symmetric in L"l. 



ii) ker(L AB ) = {cM\\c e R}. 



Proof. From the definitions of both L AB in (|1.2I) and Q AB , one has that for any u, v £ L?, 



Lab u - -^Ta^B / / (u'M'j /2 M s ; - uM l fu Blr )\{i - &) • n\d£.dn 
Mj Js+ Jm.3 

where M' B , = M B (&), Mb, = M B (f*), and u' = u(£'). 

Next notice that Ma(|£|) = M_b((^-) 1 / 2 ^) due to the energy consrevation equation satisfies 



\' A M' B . = M A M B » 

3\ 



which allows us to compute (with (,) denoting the usual inner product in R|) after the change £ — > 

- e 



(L AB u,v)=ai B I ^rj2 f (u'M^ /2 M B ;-uMy 2 M B ,)(|(e-^)-«l)^*^ 

1 /* ' 

= 9 CT ^ / ^ - ^7i7^)( u ' M a /2m s* " uMy 2 M B ,)(|(£ - 6) • n|)dnd£*d£ 
z J Ma ' 

= / m^( v ' M ^' 1/2m b* - vMY 2 M B>t )(u'M/ 1/2 M^ - ul< 2 M B „) 



- ■ n\)dnd£*dt 

which proves that Lab is self adjoint since the above expression is symmetric in u, v. Furthermore putting 
u = v we get 

{L AB u,u) = -\a\ B J ^-^(u'M^ 1/2 M^-uMY 2 M B ,) 2 (|(e-^)-n|)dr^^ 
which proves that 

(Labu, u) < 

as well as 

(L AB u,u) = uM)( 2 Mb* = u'M^ /2 M^. 

From M A M B>t = M A M' B>t and the above, one has that for all £ E W 3 

u/Mf - u'/M A /2 = 0. 

As in the argument for the Boltzmann collision invarians if for all £ we have that <j)(£) — </>'(£) = then <f> 
is constant, which with the above yields that 

u = cM| for cel. 

Thus in contrast to the usual linearized Boltzmann operator Lab has a one dimensional kernel which 
corresponds to the conservation of mass law (obtained if one integrates the original equation p. 3(1 against 
the constant 1). Notice that due to the dissipation of energy and momentum in a diffusion no further 
conservation laws could have been expected. 

□ 

Let us write L AB u as follows: 

L AB u = ^[ [ u'M A 1/2 M s ;i(£ -£*)- n\d^dn - a 2 AB [ [ uM s ,| (£-&,)• n|d£*dn 
M / J s+ Jm 3 Js+ Jm 3 

fAB.. ,,ABi 



'A 



= K AB u-is AB {\£\)u. 

We have the following lemma summarizing the above and stating some properties of the K AB — v AB 
decomposition: 
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Lemma 3.2. 

(i) It immediately follows v AB = —^-v where v is exactly as in 112.411 . so we have again 

a BB 

o< Cl <^ B (iei)<c 2 (i + iei) 

for some positive constants c\ and C2 ■ 

(ii) With the normalization RT = 1,pb = 1 the operator K^u is given by 

K AB u(£) = / k AB (Z,£,*M£,*)dt,*, where 

Jr 3 

rr 2 / m a \ 2 1 IIS*! 2 -!;! 2 ! 2 m A IS.-il 2 



which is symmetric and uniformly square integrable in each variable. 
(Hi) The operator K AB is compact L 2 (R 3 ) -> L 2 (R 3 ). 

(iv) From i, ii, and Hi we conclude that the operator Lab is a closed unbounded operator on L 2 (M?) with 
dense domain the set D = {u | (1 + |£|)u G L 2 (R 3 )} and its kernel is Ed = M^ 2 . Moreover Lab "is 
self-adjoint and non-positive, i.e.: 

(L AB u,v) = (u,L AB v) 
{L A bu, u) < 

(v) It is strictly negative on functions orthogonal to the kernel, i.e. there exists a [i > such that i/u 1 Ejj 
then 

{L AB u, u) < -m(U) u). 

The proof is analogous to the proof for the full Boltzmann linearized operator. See UJ. 



For later use we will also need the following 
Lemma 3.3. 

(i) The operator Lb a is a compact operator L'i — > L'i and is given as 



L BA }{0 = -^=Q ba (Mb,UVMa~) = [ £*)/(£* R* 

V IVI -B J 



where the kernel satisfies k BA (— — k BA (^,—^ 
(ii) The (unbounded) operator 



Lai 



'--Lba f 



1 



Q (v Myif, Mb,) 



--Q BA (M B , 



if) 



is orthogonal to all collision invariants |£| \/M~4. 

Proof, (i) The kernel k BA is obtained analogously to that in the previous lemma, (ii) To see this con- 
servative property denote by ' AB and ' BA the velocity transformations appearing in the two collisional 
integrals. Then more generally observe that 

(Q AB (f,h)(0 + Q BA (h,f)(0,m) = 



2 

a AB 



<Kf- 



lAB u'AB 



s+ 



2 

a AB 



°\b 



(L + R)(j>dnd£*d£ 
s+ 2 



/ / / (L<f> + R(j> + L<fr + R(f>)dnd£,dt. 
Jr 3 Jr. 3 Js+ 



Now notice that the transformation : £ — > , — > £ takes R to L and reversely. Moreover the 
transformation T ,j4s : £ — > £' j4B ,£* — > £,' AB takes 1/ to — L. So after performing T* on the second term, 
rp/AB Qn ^hird t erm anc j T' AB T* on the fourth term we get 



(Q AB (f,h)(o+Q BA (hj)(o,m) 
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4^ f I I m-$ AB + ^-^ AB )dnd^dt 
1 Jm. 3 Jr 3 Js+ 

From this as in the treatment of the usual collision operator it follows the above expression is if is 
a linear combination of the functions |£| 2 and the Lemma follows. □ 

Macro-Micro Decomposition for Boltzmann Diffusive equation 

Similar to the macro-micro decomposition (Po, Pi) for L, from Lemma l3~D and Lemma we can decom- 
pose L 2 into ker(LAB) © ker(LAB) ± as follows 

i = P^u + Pfu, 

P» = E D ®(E D \, (P£u = (u,E D )E D ) 

pD — i _ pD 
r l — 1 r J 

1 /2 

E D = M A ' , (ker(L AB ) = span{E D }). 



(3.1) 



3.2. The Spectrum Decomposition of Boltzmann Diffusive equation. 

Taking the Fourier Transform in the variable x of the equation (|1.3|l we obtain the equation 

8 t g = (Lab - ik^g = K AB g -vg- ik^g 

which leads us to study the spectral properties of the operator Lab — ikt; 1 . 
Similarly to the explanation in we have the following lemma. 

Lemma 3.4. For each fixed k £ R, the spectrum of Lab — ik^ 1 in the domain ReX > —vq consists of 
isolated eigenvalues with non-positive real parts. 

Next we want to consider the behavior of the spectrum for |fc| <C 1. We have the following lemma: 

Lemma 3.5. 

(i) Given kq > 0, if \k\ > k$ > 0, k G R then there exists e(«o) > such that if X is an eigenvalue of 
Lab ~ ikt, 1 then ReX < —e. 

(ii) There exist k and S such that for \k\ < k the spectrum with |A(fc)| < 5 is a single point which 
converges to the origin and it , together with its eigenvector erj(k) is an analytic function of k for \k\ < Kq. 
(Hi) The expansions of the eigenvector and eigenvalue are as follows: 

(3.2) A(fc) = —0,2k 2 + 0(fc 3 ) where 02 is a positive real number, and 

(3.3) e D (k) = E D + e' D (0)k + 0(l)|fc| 2 , where e' D (0) = iL^ l E D . 

Proof. Statements (i) and (ii) follow from standard perturbation theory and the fact that i^ B (|£|) < 
c(l + |£|), see 0. 

For each fixed k, the eigenvalue problem for the operator —ikt; 1 + Lab is 

(3.4) {-iki 1 + L AB )e D {k) = X(k)e D {k). 

Now apply the Macro-Micro decomposition (Pq , Pf) given in 13.1fl to both l|3.4|l and en(k) to get 

'eo(fc) = Poe D (k) + Pf e D {k) = ^0 + fa, 
ikPgfifa + fa) = A^ , 



(3.5) 



L AB fa + ifcPf e(i>o + fa) = Xtf) 



The last equation allows us to express 
(3.6) fa = ik[L AB - ikP? C 1 - A]-Vo- 

Which substituted into the second one in (|3.5|l with the property (E^^E^) = gives 

k*P$?[L AB - ikPi? - A]"Vo = A^o- 

1 Ii 

Taking inner product with we get the following implicit relation for X(k) 



A(A, k)=f E D {k 2 e[L AB - ikP°e - A]- 1 - X)E D d£ > = 
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We note that 



9 fe A| A)fc=0 = 0, a A A| Aife=0 = 1 
d 2 k A\ Xk=Q = J eE D L^eE D dt = -a 2 <0 



Therefore by the implicit function theorem we get 



dk 2 



dX 

dk 

k=0 



d k A 



k=0 



flgA 



= 0. 

A=0,fc=0 

= a 2 > 0. 



A=0,fe=0 



Next we consider the eigenvector e D (k) near k = 0. Normalizing e D (k) = Ejj + tl>i(k) from the expression 
(|3.6|) after differentiation we immediately obtain the desired expansion. □ 

Note that the number a 2 = — (^Ed, L^ b ^ 1 Ejj) is realized as cross species diffusion coefficient. 

4. On the Linear Solution Operators for the Boltzmann Diffusive equation 
Consider the initial value problem 
, 41 v {dtg + ^dxg = L AB g, 

^g|t=0 = gin, 

where gj„ satisfies (|1.6f) . Then, we have 

g(M) = e ( -« lfc+ ^ B)t gi„(fc). 
Now taking inverse Fourier transform as before we have for the Green's function 

<G AB {x,t) = — f e lkx+{ ~ iek+LAB)t dk. 
2tt J r 

Similarly to the Long Wave-Short Wave decomposition 12.9(1 . we also use the same decomposition for 



&AB;L{x,t) = 



B (x, t) = &AB;l{x, t) + G A B-s(x, t), 
2"7T J\k\<K 

. 8 (x,t) = -^= I e ^+(-*ek+L AB )t dk ^ 

V2n J\k\>K 

Furthermore, in the case of Gab, if ITf is the projection on the eigenvector e£>(k) corresponding to the 
eigenvalue A(fc) discussed in Lemma 13.51 we can write 

1 



"AB 



/2tt 



M-^k+L AB )t {u D + U D± )dk = G fa t) + GAB:L.X(X, t). 



I fe I < K 



Similar to 1)2. 8|) . the operator satisfies that for |fc| < Kq, 

(4.2) \\e^ lk ^ +LAB)t Ii^\\ L 2 < Ce~ t/C for some C > 0. 
Writing e D (k) = E D + e' D (0)k + 0(l)k 2 , we have that 

nf = E D <g> (E D \ + ke' D (0) <g> (Er>| + kE D <g> (e' D (0)\ + 0{l)k 2 . 
^From property (iii) of Lemma 13.51 one has that 

(4.3) et-^+^nf = e" A W* (E D ® (E D \ + fc(e' D (0) 
Here in l|4.3|l . all functions are analytic around k — 0. 



-D 



E D ®(e' D (0)\) + O(k 2 )) 
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Theorem 4.1. Recall that we denote \\\.\\\ := \\-\\L%>(L,°° 3 ) an d iPo > Pi ) the macro-micro decomposition 
given in <|3.1fl • The Boltzmann Diffusion Green function satisfies the estimates: 
In the region \x\ < 2|Ai|f, 

(x-A 2 f) 2 

(4-4) \\&ABgin(x,t)\\L2 =0(l)|||g in ||| I 6 + e -(«+H)/c ) . 



( (x-A 2 t) 2 
g (l+fl +e " (; ' ' ' ' 

( (x-A 2 i) 2 

In i/ie region \x\ > 2|Ai|£: 

||GABg m (x,i)|| L 2 = 0(l)|||g m |||e- (t+N)/C 

Proof. First, we claim the Particle- Wave decomposition in (|2.10() is also valid for l|4.1|l . since the operator 
K AB for Lab and K for L both have the properties necessary for the proofs in |13) for obtaining Lemma 
EOl Thus we have that there exist (Ff, Wf) satisfying 

T>D(~. jA i \wiDi 



(4.7) 



g(x,t) = Ff(x,t) +Wf(x,t) for all x e K,t > 0, 
Pf(x,0)=g in (x), 

Il^f 0M)llz, 2 < Cje-^+V/Cj for all x e R,t > 0, 



| (ft + ed x - Lab)*? || fl i (iS ) < Cj-e-V^i, 
l^i^llif^L 2 ) — Oj f° r a ^ ^ — (uniformly bounded in time). 

Here, we can choose j = 3. 

^From property (i) of Lemma 13.51 one can have that 

(4.8) ||Gi fl;S ginlk| (i |) < 0(l)|||g m |||e- t/C for some C> 0. 

^From property (iii) of Lemma 13.51 one can have that 

(4-9) \\G AB . L g in \\m x (q) < 0(1)11 ISmlll- 

By combining Q4.7[) and (|4.9p . one has that 

(4.10) HGxBiSgin -Pf ||fl|(Lj) = ||G^ B;i gi„ -Wf |U| (i |) = 0(l)|||g, n ||| for all i > 0. 
^From 14.71) . I4.8|l . (|4.1U|) . and Sobolev's inequality, one has that 

(4.11) \\<&AB;Sg m \\ L ~(Ll) ~ °( 1 ) e "* /C fOT SOme C > °- 

Then, by the spectral property (|4.3|) and i|3.2|) . one can apply a complex contour integral similar to that 
in ^21 to obtain that 

(■-A 2 t) 2 

(4.12) \\&AB;Lgin{x,t)\\L* < 0{l f ^—L for \x\ < 2jA X jt. 

5 V / ( 1 +*) 

Hence, 1 )4.110 and 14.12|1 conclude l|4.4|l . The estimates 1)4.511 and 1)4. 6[) will follow by taking expansion of 
higher order terms in k in (|4.3I) into account. 

In the region |x| > 2 1 Ai 1 1, one can apply the weighted energy estimates in |13j to the wave component 
W3 (x,t). Then, one can show that ||Wf (x, t)\\ L 2 decays to zero exponentially in x and t in the region 
\x\ > 3|Ai|t/2. Then, this and (|4.7|) will conclude this theorem for \x\ > 2\X 1 \t. 

□ 
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Remark 4.2. Since the decomposition (Pj^Wj 3 ) satisfies (|4.7H one can show that for large times it is 
equivalent to the decomposition (<G^ B . L , G AB .s) due to (|4.11() . We call this procedure "separation of 
scales". ■ 

5. On Cross Species Mass Diffusion and Sound Wave Interaction 
Consider the initial value problem l|1.5fl with initial data given in l|1.6fl : 

\dtg + ^d x g = L AB g, 
(5-1) \d t h + ed x h = Lh + LsAg, 

(g, h)| t=0 = (gi„, Kn)- 

The solution of the first equation is G AB gin and so satisfies the estimate provided by Theorem 14.11 Next, 
we continue to study the behavior of \\{x,i) just treating Ls A g as a source term to the equation for h. 
However, the function g and the operator dt + ^dx — L have a resonance around the trajectory x = = 0. 
This makes this system very interesting. 
1) Estimate in the Finite Mach Region x < \2Xit\ 

By Duhamel's principle one can represent the solution \n(x,t) as follows 

h(x,t) = / G(x - y,t)h in (y)dy 



(x-y,t- s)L BA G AB (y - z, s)g in (z)dzdyds. 



Due to (m and the above, we have the following "microscopic cancellation" 

1 



(5.2) 



£baJ = 



Q BA (M B ,VlvUPfJ) = L BA P»] 



D 



With this microscopic cancellation, (|4.5|) . and Theorem 12. 71 one has that 



(x~y,t~ s)L B A&AB(y - z, s)g in (z)dzdyds 



t-t2 JR 2 



/ 1 / G(x-y, t - sjLBAPi&ABiy - z,s)g in (z)dzdyds 
Jt-t? Jr 2 



<0(i)£/ > 



V 1 + (* - s ) 



+ e -(|x-J/| + (i-s))/C 



+ e -(\y\+s)/c dyds 



3 / (x-Xjt)* 

e 



+ e -(\x\+t)/C 



Furthermore with the microscopic cancellation l|5.2|) . one can rewrite Ls^Pfg in (|5.1|l through dtg + 
^dxg - L A Bg = 0: 



(5.3) 

This gives 



l ba p? g = l ba l- ab p? (d t + ea x ) 6 



(5.4) 



/ / G(x-y,t- s)L B a < Gab{v - z,s)g in {z)dzdyds 
Jo Jr 2 

i 

/ / G(x~y,t- s)L BA L' AB P°{d s +S > 1 dy)GAB(y- z,s)g in (z)dzdyds 
Jo jr 2 
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< 



+ 



I d s G(x -y,t- s)L BA L AB P®G AB (y - z, s)g in (z)dzdyds 
d y G(x -y,t- s)L B aL A b P i £ 1(G U-b(2/ - z, s)g in {z)dzdyds 



+ 



G(x -y,t- s)L BA L AB P^<G AB (y - z, s)g in (z)dzdy 



a=t* 




+ e -{\x-y\ + {t-a))/C 



+ e -(ifi+«)/c dyds 



-{\x\+t)/C 



-i\x\+t)/C 



where we noted that LBA^^sPf is a bounded operator on L? since L AB Pf is bounded, and -Lba is 
compact. Also as is clear by the proof of Theorem 12 . 71 we can justify the differentiation of G since in the 
finite Mach region G equals Gl up to a term with exponential decay in time and space. Now in the finite 
Mach region Gz, is a smooth function and the estimates on the derivatives follow from the explicit formula 
used in the proof of Theorem 12. 71 

Denote the remaining integral by q(x,t): 



(5.5) q(x,i)gi. 



t2 
1 

t — t2 „OQ 



G(x -y,t- s)L BA P?<& AB (y - z, s)g in (z)dzdyds 

^x e (- lk e+L)(t- S ) LBAP D e (- l ke+L AB )s gm ^ dkds 

For this function q(x,t), we consider the two decompositions: 
(5.6) 



t— *2 r oo 



e lkx q(k,t)dkds 



q{x,t) = qL(x,t) + qs(x,t), Long Wave-Short Wave Decomposition, 
q(x,t) = q&{x, t) + qw(x, t), Particle- Wave Decomposition, 



where 



(5.7) 



rt-sfi 


1 e ikx 


JVi J 


|fc|<K 


q(x,t) - 


qi(jc,*), 


rt-Vt 








q(x,t) - 


- q<?>(aj,i) 



J t V P ^ - s ) * ^BAPf Wf (s) + Wj(t - s) * L BA pv Pf (s) + Wj(t - s) * L BA Pf Wf da 



where 



-j, „ 3 emu i??, are given in (|2.12|) and (|4.7|l and the convolutions are in the space variable. 



and 



The consequence of these two decompositions is that for some C > 0, 

||q^>(x, t)|| £ 2 < Ce~(\ x \ +t ^ c as a convolution of decaying exponentials in time and space, 
|q^ll// 3 (L 2 ) 5; C(l + 1) as a convolution with a regular function bounded in time, 



(5.8) 



qi|li 2 (L2) < C(l +£) due to the spectral properties , 
||qs||z,|(z, 2 ) < Ce~( 1+t ^ c due to the spectral properties. 
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Combining the representations in 15. 61 and using 

f ||qs - q^||i?3(z,|) = \\qw ~ <\l\\hklI) < C(l + t), 

\hs - q^lUj(i|) < CeT l l c 
Hence, this and Sobolev's inequality give 

||qs||ig°(i,2) < 0{l)e- (1+t)/c for some C > 0. 
which in the finite Mach region x < \2Xit\ gives the desired 

(5.9) \\q s (x, t)\\ L 2 < 0(l)e~( t+ W>/ c for some C> 0. 

For < k we consider a spectral decomposition of e^ lk ^ 1+L ' ,< - t ^ s ' > LBAPi > ^~ lk ^ 1+LAB ' >s m terms of the 
spectral decompositions of the operators —ik^ 1 + L in (|2.7|) and —ik^ 1 + L A b in <|4.3I) as follows 

e (- l fcC 1 + L)(t- s ) jLs ^pZ5 e (-ifee 1 +L AB ) ;s 
3 

= JViW(t-»)+A(*)« ej .(fc) <g> (ej-(fc)|iflAPf e D (A) ® (e2,(*)| 

3=1 
3 

+ ^e^ (fc)(t - s) e J (fc)® (e J (fc)|i BA Pfe ( - 4fc « 1+L - 4B)s nf i 

+ e A(/=) Se (-«£ 1 +L)(t- S ) n ± iBA pD eD ( fc) (ec(fc) | 

+ e (- 1 ^ 1 H)( i - S )ni ijM po e (-*e 1 +i-) S nf i 

^From (E3 and JS^Jl, it follows 

(5.10) e,(fc) ® (ej(A)| L BA Pf eii(fc) ® M*)| = fcE^ ® (E,-| L flA Pf e' D (0) ® (E D | + 0(l)fc 2 

= -tfcEj ® (Ej | L SA L^Pf ?E D ®{E D \ + 0{l)k 2 . 

Remark 5.1. The components e ^( fe )(*- s )+Mk)s ei (fc) $ ( ej -(jfe)| Ls A Pf e c (fc) ® (e£>(fc)| for j = 1,3 are the 
Fourier transformation of the interaction of sound wave and cross species mass diffusion. 

The component e cr 2{k){t~s)+\(k)s e2 (/j) g, (e2(fc)| LBA^i^oik) <g> (er>(k)\ is regarded as the Fourier trans- 
formation of the "resonance of cross species mass diffusion" . ■ 

Now we will use Lemma [3.31 Cii). Let $ deonte any linear combination of l,£j, |£| 2 . We have that the 
operator 



Lab-) 



--LbaJ 



i A 

|2„ 



is orthogonal to all collision invariants ylVU, £,iVM~A, |£| 2 \/MX With this applied to J = L ab ^\Ed (where 
we recall that Ed = Ma) it follows 

J (l ab (L ab ^E d ) + ^^Lba(L ab CiE d )^ <S>VM A d£, = 



so that 



or so in the case of our interest 



E 2 ®(E 2 \L B AL AB P?e^D = 



(5.11) 



= E 2 / (- V3/2xo + Xi)L B AL- AB eE D d£, = -E s 



c « 
since the integrand is odd. This is an expression of the conservation of mass, energy and momentum in 
the entire system. 



18 



ALEXANDER SOTIROV AND SHIH-HSIEN YU 



Remark 5.2. We could have proved the above also using the reflection preserving properties of the kernels 
and the collision frequency which is implied by Lemmas 13.21 and 13 . 31 for both the operators Lab and Lb a 
Considering the transformation £ — > — £ due to the properties 

v AB {-0 = v AB {i), and k AB {-^) = k AB {i-^) 

if f(~0 = -f(0 we get 



L A±s f{-o = j k^i-t^fajdz. ~ k-o/c-o - j k AB (z,-z*)f(z*)dz*+mm = -^/(o. 

and similarly for Lb a- And similarly for an even function. From this it follows that the operator LbaL A b 
preserves the oddness of a function in £. ■ 

From (|5.10|) and 1|5.11[) 

3 

(5.12) Y.^^'^^^® (ej(k)\L BA P?e D (k)® (e D (k)\ 
i=i 

= e ^ (fc)( *- s)+A(fc)s {-ikEj <g> (Ejl L SA Pf E D ® (E_d| + 0(l)fc 2 } + e ^( fc )( t - s )+ A ( fe ) s O(l)fc 2 . 

ie{i,3} 

Now from l|5.12|l . (|2.8|) . and 14.211 we can write after performing the time integration 



(5.13) q(x,t)= f i ^- lke+L)i - t - s) L B Ae { ~ lke+LAB)s ds 



_ i U e \{k){t-t2)+a j {k)t2 _ <T J -(fc)(t-t3)+A(fe)t3 \ .2/ (A(fe)-o- 2 (fc))(t-t2) _ j\ 

H " ^ TTT^ '-^(k)+e ,J2 ^ t V ^_ ^ ^ 2 (fc)+Jf(M), 



ie{i,3} 



A(fc) - <rj{k) 



(a 2 (k)-X(k)) 



where ffj, j = 1, 2, 3, are analytic functions around k — and independent of t. 
The function J^(k,t) is an analytic function given by 



t-*2 3 



(5.14) Jf(k,t)= I ^c'iWMe^)® (e J (fc)|L Bvl Pf e (- <fc « 1 + I '^) s n^g in rf S 



j'=i 

i 

1 

/ -' 



= (_^i +L )( t _ s ) _L 



n^L BA Pf e A(fc)s e D (fc) ® (e D (fc)| g m ds 



e (^e+LKt-s) n ± LBAp D e{ -^+L)s n D± gmds _ 



Due to the spectral properties l|2.8(l and l|4.2|l . one has that for \k\ < kq 

\\Jf(k,t)\\ L 2 < 0(l)e~ tl/2/c for some C > 0. 

This yields that 



(5.15) 



f e ikx Jf(k, t)dk 



< 0(l)e 



-I* I c 



From the analytic property of the eigenvalues X(k) and Oj(k) and the expansion properties in H3.2(l and 
()2.5|l , one has that for j = 1 , 3 

(5 ' 16) A (fc )-t j(fc ) ^ (fc) = ^WTW^ (fc) 

is an analytic function of A; for k < kq which is independent of t. ^From l|5.16[) . by shifting the contour of 
integration using the analyticity as done in J3|> one has that for |x| < 2|Ai|t 
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(5.17) ^ 



. j te \(k)(t-t?)+a j (k)t? 

\k\< K0 J \(k)-aj(k) 

f ie u j (k){t-ti)+\(k)t'S 



Akx 



|fc|<K 



\{h)-aj{k) 



fj(k)dk 
7 3 (k)dk 



p C(t + 1) 

< Oil) , for j = 1,3 

(x-A.t) 2 



e c(*+i) 
< 0(1)—==- for j = l,3 

V^+t) 



We just mention that the new contour of integration is of the form 51=riUr2Ur3 where Ti = (— kq/2, — kq/2^- 
ic), T2 = (—kq/2 + ic, no/2 + ic), T3 = (kq/2 + ic, kq/2). Here the constant c < kq is appropriately chosen 
as in [T3] . 

Next we can write 

(a 2 (fc) - A(fc)) 

r 1 1 

= e (-ia 2 fe 2 + 0(fe 3 ))t fc 2^ _ t ^ e (-i a2 k 2 +0(k 3 ))t / e (A(fc)-<r 2 (fc))(t-t5)u du ^ 2 ^ 

.' i 

= e (-i a2 ^ +0 (fc 3 ))t fc 2 (t _ ^)^(fc, t) 

where (?<j;(k,t) is clearly an analytic function of k and can be easily seen to be bounded in t along the 
contour £1 by a constant independent of k and x provided k is sufficiently small and \x\ < |2Aif|. This 
allows one to again shift the contour of integration as in ^] and obtain that for \x\ < 2|Ai|i 



(5.19) 



h 2l (A(fe)-cr 2 (fe))(t-t5) _ ,\ 

e'to e ^W« - KG , — r"7TT\ -J-0 2 (k)dk 



\k\<K 



(a 2 (k) - \(k)) 



p C(t+1) p C(t+1) 

<tO(l) 7 ___ = 0(l)- 



(l + <) 3 / 2 



(l+t)V2 



Finally, from the estimates in (|5.15|l . (j5.17|l . and H5.19[l combined with the spectral decomposition of 
q(fc) in (|5.13|) and the definition 1)5. 7fl one has that 



(5.20) 



q L (x,t)\\ L , < 0(1) j ^ e ~^== ) +^ ti/C ) for N < 2|A 1 |(t + l). 



Now from l|5.20|l . I]5.9p. (|5.3fl . and (|5.4|l we conclude the estimate in the region |x| < 2|Ai<|: 



3 



e c(i+t) 



(5.21) \\h(x,t)\\ L2 <0(l) X' +e~ t§/C 



2) Estimate in the region |x| > 2|A]t| 
We simply use Duhamel's principle, 



(5.22) \\h(x,t)\\ L >< 



G(x - y,t)h in (y)dy 

&{x -y,t- s)L BA &AB(y - z, s)g in (z)dzdyds 



e c<i+t) 



•0(1) 



(x-j,-Aj(t- s )) 2 



5 ^/fi 



J V(l + (*-a)) 



3 

+ e -(k-»l+(t--))/o I I ^ e " 



-(|j/l+s)/C 



dyds 
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ALEXANDER SOTIROV AND SHIH-HSIEN YU 



< (i) y e 2C(1+<) + e -(M+t)/[2C) for | s | > 2 | A | (< + iy 

Finally, (JOJl, (|OT|) . f5^|) conclude Theorem O 
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